Frontiers in Computational Neuroscience

Cerebellar neural circuit dynamics rely on a rich repertoire of synaptic and excitability mechanisms, which are thought to determine network computation in physiological and pathological conditions. In this work, we develop and validate a biologically-grounded spiking neural network of the cerebellar cortex, embedding key mechanisms of cellular excitability and synaptic transmission, and assess their impact on signal processing. Neuronal input-output functions, short-term synaptic plasticity, receptor-specific kinetics, and NMDA channel voltage-dependent gating were calibrated against detailed multicompartmental models through automatic tuning procedures. Incorporating these realistic biological properties allowed the network model to simulate key features observed in recordings from acute cerebellar slices. The neuronal discharge and local field potentials elicited by mossy fiber stimulation faithfully reproduced the natural patterns with millisecond precision. Then, selective receptor switch-off revealed the contribution of NMDA, GABA, and AMPA receptors to the frequency-dependent input-output function of the granular layer and Purkinje cells, linking previous empirical findings to specific synaptic mechanisms. This model combines high computational performance with biological realism and offers a computationally efficient framework to investigate neurophysiological phenomena and the neural correlates of behavior in large-scale long-lasting simulations, such as those needed to address the neural underpinnings of learning and of cerebellar pathologies.

Categorisation is often treated as a form of compression: a high-dimensional stimulus space is reduced to a smaller set of behaviourally or cognitively useful classes. However, compression alone does not determine whether a category map is useful. The present manuscript develops an information-theoretic framework for evaluating categorisation in terms of both category complexity and target-relevant information preservation. Across a set of synthetic demonstrations, alternative category maps over the same stimulus space are shown to preserve different target variables, including identity, action, nuisance, and hierarchical category structure. The framework is then extended to learned visual representations by analysing layer-derived category maps from a pretrained ResNet-50 network applied to CIFAR-10 images. Two scenarios are compared: a clean-only object run and a pooled nuisance run containing clean, blurred, pixelated, and noise-perturbed images. The results show that category maps can have substantial entropy while preserving information about a variable that is not aligned with the specified target, and that the value of a categorisation depends on the target variable to be preserved. The manuscript argues that categorisation should therefore be evaluated not only by compression or separability, but by the information retained about a specified cognitive, behavioural, or computational target.

Equilibrium propagation (EP) is a biologically plausible alternative to backpropagation that has demonstrated competitive performance across a range of machine learning tasks. Recent work has extended EP to spiking neural networks (SNNs), leveraging leaky integrate-and-fire (LIF) neurons and spike-based plasticity rules to improve biological realism while maintaining strong performance. In this work, we propose an EP-based SNN framework that combines LIF neural dynamics with a predictive learning rule, replacing conventional spike-timing-dependent plasticity (STDP) with a learning rule more directly aligned with predictive coding principles. We evaluate the proposed model on multiple image classification benchmarks, including MNIST, KMNIST, and Fashion-MNIST, and compare its performance with a BP-trained LIF SNN baseline. Our results show that the proposed EP-based LIF model (EP+LIF) achieves competitive accuracy across datasets, with performance approaching that of the BP-trained counterpart (BP+LIF) while relying on a biologically motivated local learning rule. In addition, analysis of hidden-layer spiking activity reveals that EP+LIF produces more persistent hidden-state activity, whereas BP+LIF yields sparser spiking representations. These results demonstrate that predictive learning can support effective EP-based training in LIF spiking networks, while also highlighting differences in activity patterns that motivate future work on activity regulation and sparse spiking dynamics.

The awake brain is known to display asynchronous (AS) states during periods of attention and arousal, but the responsiveness properties of such states remain unclear. Here, we investigate this question using computational models of spiking networks of excitatory and inhibitory neurons, mimicking recurrently-connected networks in layer 2/3 of the cerebral cortex. The networks can generate a continuum of AS states, but with different responsiveness characteristics. By using a mean-field model to infer the dynamic properties of the system, we find that there are two families of AS states, which we call "underdamped" (UD) and "overdamped" (OD). Responsiveness is maximised at the transition between OD and UD states, which identifies a "working point" that may present advantageous computational properties.

The learning dynamics of biological brains and artificial neural networks are of interest to both neuroscience and machine learning. A key difference between them is that neural networks are often trained from a randomly initialized state whereas each brain is the product of generations of evolutionary optimization, yielding innate structures that enable few-shot learning and inbuilt reflexes. Artificial neural networks, by contrast, require non-ethological quantities of training data to attain comparable performance. To investigate the effect of evolutionary optimization on the learning dynamics of neural networks, we combined algorithms simulating natural selection and online learning to produce a method for evolutionarily conditioning artificial neural networks, and applied it to both reinforcement and supervised learning contexts. We found the evolutionary conditioning algorithm, by itself, performs comparably to an unoptimized baseline. However, evolutionarily conditioned networks show signs of unique and latent learning dynamics, and can be rapidly fine-tuned to optimal performance. These results suggest evolution constitutes an inductive bias that tunes neural systems to enable rapid learning.

Improving long-term memory in artificial neural networks remains an open challenge. To address this, we developed a novel brain-inspired framework for memory prioritization based on the principle of emotional valence. Our framework includes: (i) a valence-weighted cross-entropy loss that scales the learning signal by the valence magnitude, analogous to neuromodulation; (ii) an amygdala-inspired module that learns high-dimensional valence embeddings; and (iii) a hippocampus-inspired module that integrates valence embeddings into the attention mechanism to modulate information retrieval. We demonstrated the generalization of our framework across spatial, episodic, and language-based memory tasks, consistently improving memory prioritization and long-term retention of high-salience information. In addition to improving long-term memory, we also showed that our framework can help mitigate the "lost-in-the-middle" problem in language modeling. More generally, this research provides further evidence of the potential of brain-inspired algorithms to advance the field of machine learning.

The routing of information across the brains structural network is central to its wide range of functional capabilities. However, the mechanisms underlying information routing in complex brain networks, particularly between regions that do not share a direct anatomical connection, remain poorly understood. Neural mass models (NMMs), a computational modelling framework capable of capturing complex neural dynamics across scales, can potentially be used to study the dynamical and network bases of these vital polysynaptic routing processes. In this study, we investigate polysynaptic signalling in three widely used NMMs, obeying Ornstein-Uhlenbeck, Stuart-Landau, and Jansen-Rit dynamics, by tracking the propagation of a discrete, focal, high-amplitude perturbation across the underlying network. We find that polysynaptic propagation emerges in all tested NMMs when configured within dynamical regimes that effectively enhance the persistence of perturbations. We also find distinct parameter domains that maximise signal propagation to directly connected regions and to those separated from the source by at least two hops. Finally, we benchmark in silico stimulus propagation in the brain network against an empirical dataset of direct electrical stimulation trials, to explore the relative capabilities of the NMMs in capturing signal propagation to connected versus unconnected regions. This analysis highlights the significance of dynamical repertoire in capturing stimulus propagation outcomes. Overall, this study provides insights into how dynamical and network features shape signal propagation over complex brain networks.

Adolescence is an important developmental period during which there are diverse changes in the brain and behavior. Goal-directed behaviors and the component processes underlying those behaviors improve during adolescence, including working memory, response inhibition, and reinforcement learning. At the same time there is substantial pruning of excitatory connections in prefrontal cortex and ongoing myelination of axons. However, psychiatric disorders also become increasingly prevalent in late adolescence and early adulthood. In this study, we develop computational models that suggest a hypothesis for how the ongoing changes in the brain can give rise to the increased prevalence of psychiatric disorders. We show that both myelination and pruning during adolescence lead to attractor landscapes in which strongly encoded memories, driven by three-factor learning rules that modulate Hebbian plasticity, come to dominate the landscape of brain activity, at the expense of weakly encoded memories. Pruning and myelination lead to large, strong attractors which, if they are related to aversive emotions, can drive intrusive thoughts and compulsions in obsessive compulsive disorder, rumination in depression, and aversive memories in post-traumatic stress disorder. The link between pruning, myelination and the emergence of dominant attractors for emotionally salient memories is well supported by the models. The way these effects map onto forebrain circuits requires more work.

The phenomenon of splitting was originally observed in hamsters which, after prolonged exposure to constant light, exhibit two rest/wake cycles within a subjective day. Splitting is a consequence of the left and right suprachiasmatic nuclei (SCN) falling out of synchrony. While it is known that split activity is characterized by an antiphase relationship between the left and right SCN and between the core and shell within each hemisphere, the role of the commissural projections that connect the right and left SCN is not known. In the present study, we investigate the impact of the inter-hemispheric connections on the split and unsplit dynamics of a computational model of the bilateral SCN. Our model has 4 nodes corresponding to each right and left core and shell. We simulated our bilateral model under different lighting conditions and measured its period and the phase relationships among the 4 nodes. To further characterize the dynamics of the system, we performed a bifurcation analysis. We found that the bilateral model automatically splits unless entrained by bright light/dark cycles, or unless it has excitatory inter-hemispheric connections. This suggests that excitatory cross-connections may be important for freerunning behavior. We found that constant light of varying intensities transitions the model between split and unsplit activity only in very limited conditions, but the strength and polarity of the contralateral connections play a much greater role in this dynamical transition. These findings suggest that splitting may involve plasticity of the inter-hemispheric connections of the SCN.

In previous studies by the author on binocular vision with the asymmetric eye (AE), which models a healthy human eye with misaligned optical components, the results were primarily presented in the Rodrigues vector (RV) framework and supported by simulations and 3D visualizations in GeoGebras dynamic geometry environment. In this paper, the novel geometric kinematics of the human eye, that is, the eye with misaligned optics, and simplified assumptions about the eye rotations (the eyes translational movements are disregarded), are developed within the framework of rigid-body rotations. The originality of the analysis lies in a precise geometric decomposition of a full rotation of the eyes posture into a torsion-free rotation (the geodesic part) and a torsional rotation (the non-geodesic extension of the geodesic part). This decomposition is extended to the corresponding decomposition of the angular velocity. A novel derivation of the eyes angular velocity from the RV formulation of the eye kinematics is proposed.

Growing interest to describe the electrical behavior of glial cells, mainly astrocytes, in intact brain tissue poses more and more challenges to commonly accepted belief they only respond in a linear manner in uptake of the excess of extracellular potassium and maintenance of their network equipotentiality. Their highly conductive mutual interconnections via gap junction (GJ) connections introduce yet another class of nonlinear elements. As more studies report nonlinearities in membrane voltage Vm dependence of both, the membrane and junctional conductances, the need to formulate minimal dynamical models of their transient behavior is getting more acute. Since ODE models of coupled cells, even in simplest 1-d arrays, require simplified descriptions and small set of parameters, rare quantitative studies on glia makes the task even more difficult. This study attempts to qualify a self-coupled cell, or a glial cell coupled to fixed voltage as useful system for detecting the nature of instabilities and transitions coming from coupling. In a novel biophysical model of coupled astrocyte, we introduce nonlinear kinetics of deactivation for large junctional voltages for the first time. We found that N-shaped nonlinearities and corresponding fold structure in the vector field of isolated cell serves as a baseline on top of which coupling nonlinearities enrich the bifurcation picture. Numerical simulations of 1-d array of coupled astrocytes show that coupling increases the propensity of astrocytic Vm to bistability and front propagation. We believe that presented illustrations of possible effects of coupling nonlinearities will motivate neurobiologists to further explore their impact in disease. Significance statementTransient changes in membrane voltage of glial cells may produce significant transient voltage difference between directly coupled cells. Nonlinear steady-state conductance of their interconnection elements, the gap junctions, introduce nonlinear current profiles which are very difficult to measure and quantitate using the available methods due to marked permeability of the junctions and leakiness of glial membrane in general. We propose a minimal model of glial membrane extended with a self-coupled feedback loop, which under realistic simplifying assumptions could serve for qualitative analysis of the impact of coupling, on the stability of resting membrane voltage. Neuronal cells of the brain and spinal cord cannot exist and function without supportive and neuromodulatory functions of the diverse population of glial cells. This applies to virtually all physiological processes on cell level - from cell development, metabolic support, membrane signaling, slow molecular signal transduction, ion homeostasis, neurovascular coupling, myelination, to mention only a few, manifest neuro-glial interaction. Even though all glial cell types are interconnected, the most abundant ones, the astrocytes are massively interconnected by gap junctions to form ordered networks. Electrically, astrocytic networks display membrane voltage equipotentiality, which is considered system-wide resting state for given neuro-glial circuit or unit. With molecular and cellular substrates of glial connectivity being slowly elucidated, network science and dynamical modeling are slowly "invading" that area with many important issues left open. In this study using classical dynamical systems approaches we give indications how nonlinear intercellular coupling between astrocytes affects physiological resting state and its instabilities compared to isolated, uncoupled cell. We strongly believe the suggested minimal model could fill the gap in ODE modeling of neuro-glial circuits, within broadest scope of hypothesis-driven research in cell-level neuroscience.

Rhythm is a key building block of human music, speech and numerous other human activities. Understanding the computational substrates of rhythm perception requires models that bridge algorithmic function with biological implementation. We propose a physiologically grounded spiking neural network (SNN) framework to investigate the emergent representation and interpretation of auditory rhythms. Utilizing a recurrent SNN architecture trained on an auditory entrainment task, we characterize the networks latent dynamics through the analysis of firing rates and membrane potential fluctuations. Our results demonstrate that simulated neural populations exhibit phase-locking to the stimulus beat, with endogenous oscillations driven by rhythmic input. We further show that anticipatory dynamics--characterized by pre-stimulus depolarization--emerge naturally from the networks synaptic plasticity and temporal integration properties, rather than from explicitly defined oscillators. By treating network layers as functional analogs of cortical populations, this framework allows for the application of spectral and information-theoretic analyses typical of empirical electrophysiology. More in general, this approach establishes SNNs as robust exploratory tools for uncovering how predictive coding and rhythmic entrainment arise from the inherent constraints of biological neural computation.

We find that multi-site temporal control of optogenetic photostimulation in peripheral nerves can enhance firing rates by overcoming the intrinsic limitation of opsin photophysics. The benefits of multi-site optogenetic stimulation were demonstrated with three approaches: (1) in silico modeling, (2) ex vivo in the sciatic nerve, and (3) in vivo in the vagus nerve. An in silico model of multi-site optogenetic stimulation was developed in two Hodgkin and Huxley type neuron models, that supported our hypothesis. The ex vivo sciatic nerve showed an increase in firing frequency that is physiologically relevant for functional control. The technique was then applied in vivo for optogenetic vagus nerve stimulation resulting in significant changes in heart rate compared with standard methods of single-site stimulation. Improving the control of optogenetically induced neural firing will have broad impacts for future developments in optical nerve interfaces and brain-machine interfaces.

Metabolic networks are typically viewed as homeostatic systems that stabilize flux, energy charge, redox balance, and metabolite availability under perturbation. However, it remains unclear whether the same feedback architectures that support metabolic robustness can also generate learning-like, experience-dependent adaptation. Here, we develop a coarse-grained dynamical model of mammalian energy metabolism to test whether prior perturbation can improve future metabolic responses. The model represents core glucose, glutamine, fatty acid, and oxidative phosphorylation pathways as coupled ordinary differential equations with Michaelis-Menten-type fluxes, product-inhibition feedback, adaptive enzyme-capacity regulation, and explicit ATP costs for enzyme adjustment. Rather than aiming to reproduce quantitative fluxes for a specific cell type, the framework is designed to expose how metabolic feedback, regulatory cost, repeated perturbation, and environmental variability interact. We use this model to ask whether adaptive enzyme regulation enables improved recovery after repeated challenges, whether such effects depend on energetic control costs, and whether environmental variability broadens or constrains the set of reachable adaptive states. This approach provides a tractable way to investigate how homeostatic metabolic regulation may give rise to experience-dependent metabolic plasticity.

Humans frequently make decisions in complex, high-dimensional environments, where identifying task-relevant information is critical for rapid behavior optimization. Humans outperform standard reinforcement learning agents in navigating such complexity, yet the cognitive strategies of humans remain unclear. To address this, we developed a novel multi-dimensional learning task in which only a subset of dimensions is reward-related. Crucially, unlike prior studies, subjects are uninformed of the true task dimensionality and have to identify them through exploration. This design closely mimics the ambiguity in real-world tasks. Our results have identified two stereotyped choice patterns that reveal "dimension-guided" strategies in exploration and exploitation. Cross-subject analyses suggest that dimension-guided exploration may promote the efficiency of reward-based learning. These findings indicate that humans leverage task dimensionality to guide exploration, and provide inspiration for improving exploration efficiency in AI agents.

A growing number of studies indicate the possible involvement of astrocytes in triggering or modulating neurovascular coupling (NVC), i.e. the local dilation of blood vessels in the brain in response to neuronal activity. Astrocytes possess specialized subcellular compartments, named endfeet, that surround arterioles and capillaries, ideally positioned to mediate NVC. Various vasodilators have been shown to contribute to NVC, such as epoxyeicosatrienoic acid (EET), nitric oxide (NO), or prostaglandin E2 (PGE2), but the precise mechanisms underlying NVC and their variability remain to be fully elucidated. In particular, the involvement of astrocytes in this process is controversial. Recent translatome and proteomics data reveal that astrocytes and in particular endfeet are enriched in the proteins of the PGE2 pathway. However, how the latter could contribute to NVC remains to be characterized. Here, we develop a computational model of astrocyte-mediated NVC that recapitulates these findings and describes Ca2+ and PGE2 signaling in astrocytes, NO release by neurons, and arteriole diameter dynamics using ordinary differential equations. The model successfully reproduces the dynamics of arteriole diameter change during hyperemia from in vivo neocortical recordings in awake mice. Our simulations suggest that the astrocyte PGE2 pathway could be responsible for the late response of NVC at the arteriolar level. We further observe that PIP2-derived diacylglycerol plays a major role in driving arteriole diameter dynamics in our model, while phosphatidic acid-derived diacylglycerol, which is calcium-dependent, mainly acts as an amplifier of this response. Finally, a spatial implementation of the model using a simplified astrocyte geometry suggests that NVC is more efficient when synaptic stimulation occurs at the endfoot level rather than at other astrocytic compartments. Overall, this computational study suggests a partial role for astrocyte-mediated PGE2 release in NVC and points to astrocyte perivascular processes as sub-compartments that are ideally positioned and equipped to mediate NVC. Author summaryIn the brain, the local blood flow is regulated to meet neuronal energy demand by modulating the dilation of neighboring blood vessels. The mechanisms driving this process, known as neurovascular coupling (NVC), remain debated and are likely to differ depending on the physiological context. Recent evidence points to astrocytes, a cell type possessing specialized protrusions called "endfeet", that envelop the entire brain vascular tree. Contacts between synapses and endfeet have recently been reported, positioning the latter as ideal mediators of NVC. Here, we developed a computational model that simulates the signaling between neurons, astrocytes, and blood vessels. Our model successfully reproduces experimental recordings of blood vessels dilation in the brains of awake mice. Our simulations suggest that a specific signaling pathway in astrocytes, involving a molecule called prostaglandin E2, is a key driver of the late phase of NVC, occurring a few seconds after neuronal activity. Furthermore, our model indicates that the location of the stimulated synapses matters: signals sent to the astrocyte endfeet are particularly effective at controlling blood flow. This work helps clarify the active role of astrocytes in brain blood flow regulation, a process critical for healthy brain function.

The acquisition of sensorimotor skills critically depends on basal ganglia (BG)-thalamo-cortical circuits. Prevailing theories propose that the BG optimize motor output through reinforcement learning (RL), using internal performance evaluations to approximate stochastic gradient ascent. However, this framework struggles in non-convex performance landscapes, where local optima hinder efficient learning. Songbirds provide a compelling biological example of robust sensorimotor learning, mastering complex vocalizations through trial-and-error within a specialized BG-thalamo-cortical architecture. Here, we present a computational model constrained by the anatomy, physiology, and developmental trajectory of the zebra finch song system. The model combines a BG-driven RL pathway with a parallel cortical motor pathway that progressively consolidates successful motor patterns via Hebbian plasticity. In addition, we incorporate synaptic volatility within the BG pathway, introducing structured variability across learning. Through simulations of vocal learning using both a biophysical syrinx model and synthetic performance landscapes, we demonstrate that this dual-pathway architecture reliably converges to global optima and outperforms standard and noise-annealed RL approaches. The model reproduces key experimental features of song learning, including non-monotonic learning trajectories, a gradual reduction in motor variability, and the developmental transfer of motor control from subcortical to cortical circuits. Mechanistically, delayed maturation of the cortical pathway provides an implicit regulation of the exploration-exploitation trade-off, while synaptic volatility enables escape from local optima. These results highlight the importance of neural circuit architecture and dynamics in efficient learning, and suggest biologically inspired design principles for improving the robustness and sample efficiency of artificial RL systems in complex sensorimotor domains.

Navigation in complex environments relies on internal spatial representations that guide action. While the brain employs a diverse repertoire of spatial tuning cells--including grid, place, and head-direction cells--a normative theory linking these static neural codes to the dynamic process of navigation remains elusive. In this work, we propose a Unified Laplacian Framework derived from first principles of representational smoothness and efficiency. We first demonstrate that diverse spatial codes emerge naturally as spectral decompositions of the Laplace operator. Crucially, bridging the gap from representation to action, we derive a biologically plausible navigation policy based on the Greens function potential. We show that this potential encodes the environments intrinsic geometry to enable simple, trap-free gradient ascent, achieving improved sample efficiency and generalization in goal-reaching tasks. Furthermore, we demonstrate that these spectral representations can be learned directly from high-dimensional visual inputs, confirming their plausibility in realistic environments. Our results suggest that the "cognitive map" can be viewed as a spectral embedding of the Laplacian, providing a rigorous foundation for spatial cognition in both biological and artificial agents.

A motoneuron pool is often regarded as a rigid controller because the largely shared synaptic input across motoneurons leads to strongly correlated activity. However, brief deviations from this correlated behavior have been observed even in some constrained tasks, raising the question of whether these results reflect limitations of the rigid view of motoneuron pool control. Here we show that they do not. We developed a biophysical model of a motoneuron pool receiving shared excitatory and inhibitory synaptic inputs that also included the motoneuron-specific effects of neuromodulation; model parameters were tuned based on large-scale motoneuron recordings in humans. Simulations showed that the intrinsic differences in how motoneurons respond to neuromodulation are both necessary and sufficient to transiently decorrelate pairs of motoneurons receiving a shared synaptic input. Crucially, such transient decorrelation is only observed when motoneurons have different sensitivity to neuromodulation, consistent with experimental observations during volitional control in humans. Our model also explains how participants can improve their ability to transiently decorrelate the activity of motoneurons innervating the same muscle by leveraging refined behavioral strategies that exploit the differential response of motoneurons to neuromodulation, rather than through physiological changes. These results identify that heterogeneous sensitivity to neuromodulation enables brief flexibility in the otherwise rigid control of motoneurons enforced by a shared synaptic input, and show how practice allows participants to exploit latent flexibility within otherwise rigid constraints.

Understanding intracranial pressure (ICP) dynamics is essential for interpreting clinical infusion tests used in the diagnosis of cerebrospinal fluid circulation disorders. However, the complex coupling between vascular pulsations, cerebrospinal fluid flow, and intracranial compliance makes quantitative interpretation of these tests challenging. Here, I present a patient specific simulation framework based on an extended electrical analog model that reproduces intracranial pressure dynamics observed during clinical infusion tests. The model integrates physiological inputs including arterial blood pressure, heart rate, respiratory rhythm, and resistance to cerebrospinal fluid outflow derived from clinical data. Built upon the classical Ursino framework, the model incorporates several modifications enabling realistic representation of physiological pulsations and infusion test conditions. The resulting system functions as a hybrid electrical-numerical simulation model representing a simplified digital electrical twin of intracranial hydrodynamics. The model was validated using data from 21 clinical infusion tests performed in patients with suspected normal pressure hydrocephalus. Simulated intracranial pressure recordings were compared with clinical measurements using regression and residual analysis. The simulations demonstrated strong agreement with measured data, with a mean correlation coefficient of r = 0.95 (95% CI 0.94 - 0.96), mean residual values within -1.71 to +1.68 mmHg, and a mean root mean square error (RMSE) of 2.07 mmHg. These results demonstrate that the proposed model accurately reproduces the dynamic behavior of intracranial pressure observed during clinical infusion tests. The framework provides a physiologically grounded computational tool for studying patient specific intracranial dynamics and may support improved interpretation of infusion test results in clinical practice.